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1.  Introduction 


Future  military  engagements  will  require  weapons  systems  exhibiting  improved  range  and 
accuracy.  One  of  the  technologies  under  investigation  at  the  U.S.  Army  Research  Laboratory 
(ARL)  to  achieve  these  goals  is  the  electrothermal-chemical  (ETC)  propulsion  concept,  shown 
schematically  in  figure  1 .  In  the  ETC  gun,  energy,  which  is  stored  either  in  batteries  or  a  rotating 
device,  is  converted  on  demand  into  an  electrically  generated  plasma  (resulting  from  the  ablation 
of  polyethylene  material  in  a  capillary)  that  is  injected  into  the  chamber  of  a  howitzer  or  gun. 

This  plasma  energy  is  used  to  ignite  the  chemical  propulsion  charge  (i.e.,  solid  propellant)  as 
well  as  to  enhance  gun  performance  by  taking  advantage  of  a  number  of  unique  plasma 
characteristics.  For  example,  a  low  density  plasma  jet  can  efficiently  ignite  charges  of  high 
loading  density,  can  control  propellant  mass  generation  rates  (7),  can  reduce  propellant  charge 
temperature  sensitivity  (i.e.,  the  variation  of  gun  performance  with  changing  ambient 
temperature  [2,  5]),  and  can  shorten  ignition  delay  (i.e.,  the  time  interval  between  firing  of  the 
igniter  and  ignition  of  the  propellant  [ 4] ).  Plasma  igniters  replace  the  chemical  igniter  and  thus 
enhance  the  safety  of  the  gun  system. 


Figure  1.  Electrothermal-chemical  (ETC)  gun  concept. 

Research  has  been  carried  out  on  the  use  of  plasmas  to  ignite  solid  propellant  (5,  6).  Since  the 
plasma  is  at  a  temperature  (-15,000  K)  that  is  higher  than  chemical  igniters  (-3,000  K),  the 
radiation  properties  of  the  plasma  have  also  been  considered.  The  high  plasma  temperature  leads 
to  radiation  effects  nearly  100  x  greater  than  that  of  chemical  igniters  (i.e.,  a  T4  effect)  (5).  Such 
radiation  could  lead  to  significantly  different  temperature  profiles  within  the  propellant,  causing 
changes  in  bum  rates.  In  addition,  plasma  has  a  much  lower  density  than  the  gases  generated  by 
a  chemical  igniter,  a  feature  that  alters  the  convective  heat  transfer  to  the  propellant  as  the 
plasma  moves  through  the  grains  as  well  as  the  velocity  and  mode  of  flame  spreading  within  a 
propellant  bed.  It  has  been  suggested  that  energy  transport  by  convection  may  be  as  important  as 
radiation  transport  in  plasma-propellant  interactions  (PPIs)  (7,  8). 
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The  effects  previously  described  can  lead  to  significant  changes  in  ballistic  behavior  and  useful 
improvements  in  gun  performance,  but  an  understanding  of  the  underlying  physical  mechanisms 
is  necessary  in  order  to  achieve  these  goals.  To  this  end,  ARL  is  supporting  a  comprehensive 
study  of  the  interaction  of  the  plasma  efflux  from  an  ETC  igniter  with  solid  propellant.  The  goal 
of  this  work  is  to  elucidate  the  relevant  physical,  mechanical,  and  chemical  mechanisms  that 
underlie  the  observed  ballistic  effects.  Various  aspects  of  the  experimental  and  modeling 
research  program  for  PPI  are  described  elsewhere  (9-20).  The  present  report  explores  progress 
in  the  symbiosis  of  experimental  work  and  a  multidimensional,  multiphase  interior  ballistics 
code,  NGEN,  that  can  explicitly  model  the  two-dimensional  chamber  of  a  gun  (including 
chambrage),  the  solid  propellant  charge,  and  the  inflow  of  plasma  into  the  gun  chamber  from  an 
ablative  polyethylene  capillary.  The  NGEN  code  predicts  propellant  ignition  and  flame 
spreading,  chamber  pressurization,  pressure  wave  travel,  projectile  movement,  and  propellant 
burnout.  The  prediction  of  these  interior  ballistics  parameters  is  important  to  the  design  of  ETC 
guns.  Code  validation  using  data  collected  in  a  25-mm  ballistic  simulator  is  discussed. 


2.  Experimental  Description  for  the  ARL  25-mm  Simulator 


The  ARL  25-mm  ETC  ballistic  simulator  experiment,  constructed  by  Chang  and  Howard 
(21-23),  is  used  to  investigate  plasma  flow  through  the  center  core  of  an  assembly  of  concentric 
solid  propellant  disks  that  are  each  separated  axially  by  a  small  gap — much  like  the  configuration 
of  an  actual  gun  charge.  Figure  2  shows  a  chamber  made  of  visually  transparent  acrylic  that 
allows  cinematography  of  plasma  flows  and  ignition  events  along  the  propellant  bed.  The 
chamber  can  withstand  pressures  up  to  ~13  MPa  (1,900  psig)  before  rupture.  Two  Kistler 
pressure  transducers  (Model  21  IB  1)  are  installed  for  monitoring  chamber  pressures  near  the 
breech  and  forward  ends  of  the  chamber.  The  chamber  is  10.15  cm  long  and  3.4  cm  inside 
diameter  (ID).  The  chamber  is  narrowed  to  2.44  cm  diameter  at  the  forward  end  and  sealed.  A 
projectile  afterbody  of  1.5  cm  length  and  2.1  cm  base  diameter  is  inserted  into  the  forward  end  of 
the  chamber.  Along  the  length  that  is  of  constant  diameter,  a  series  of  concentric  disks  are 
placed  (figure  2).  Initial  experiments  utilize  disks  made  from  the  inert  material  Masonite*  with 
dimensions  of  33  mm  outside  diameter  (OD),  8  mm  ID,  and  6.35  mm  thickness  (5.84  g  per  disk). 
Live  propellant  disks  are  made  from  graphite-free  JA2  with  dimensions  of  33  mm  OD, 

9.4  mm  ID,  and  2.29  mm  thickness  (2.64  g  per  disk).  The  chambrage  region  surrounding  the 
afterbody  of  the  projectile  is  filled  with  7-perf  gr  of  inert  material  (2.54  mm  in  length  and 
2.54  mm  in  diameter,  perforation  diameter  of  0.33  mm)  for  both  inert  and  live  experiments. 

In  two  previous  setups,  12  disks  of  inert  material  were  loaded  with  an  initial  gap  between  each 
disk  of  about  0.8  mm,  or  13  disks  (last  disk  was  5.35  mm  thick)  were  loaded  with  an  initial 


Masonite  is  a  registered  trademark  of  Masonite  International  Corporation. 
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Figure  2.  Photographs  of  the  ARL  25 -mm  simulator  loaded  with  disks  without  spacers  and  a  single 
Masonite  disk;  dimensions  shown  are  inches  (22). 

gap  between  each  disk  of  about  0.4  mm  (21-23).  In  each  case,  the  chamber  space  ahead  of 
the  last  disk  was  filled  with  small  7-perf  gr  of  inert  material  (2.54  mm  in  length  and  2.54  mm  in 
diameter);  a  small  0.8-mm  gap  existed  between  the  breech  wall  and  the  first  disk  and  a  small 
0.5-mm  radial  gap  existed  along  the  radial  wall  of  the  chamber.  These  setups  were  used  to 
study  the  performance  of  the  NGEN  code  for  the  case  of  inert  disks  (24).  In  the  latest  setup, 

30  graphite-free  JA2  disks  with  a  0.5-mm  gap  spacing  were  used.  This  current  setup  is  studied 
in  the  present  report. 

The  ablative  capillary  used  in  this  study  to  generate  the  plasma  basically  consisted  of  a  stainless 
steel  cylinder  that  housed  a  polyethylene  tube  with  an  ID  of  3  mm.  The  igniter  wire  for  the 
plasma  was  the  same  0.13-mm-diameter  nickel  wire  used  in  prior  experiments.  The  anode  was 
made  of  tungsten  and  a  stainless  steel  cathode  formed  the  plasma  jet  exit.  The  capillary  was 
contained  within  a  plasma  injector  which  was  mounted  on  the  breech  end  of  the  chamber  (see 
figure  2).  The  pulse  power  supply  was  capable  of  delivering  energies  of  up  to  nearly  3  kJ  at  a 
voltage  of  3  kV.  The  typical  pulse  length  from  the  power  supply  was  about  360  ps.  The  pulse 
was  lengthened  to  ~1  ps  for  several  experiments  by  the  addition  of  a  44-mH  inductor  to  the 
positive  lead  prior  to  attachment  to  the  capillary.  The  voltage,  current,  chamber  pressures,  and 
ignition  events  were  recorded  directly  on  a  Nicolet  Integra  20  digital  oscilloscope.  With  the 
measured  voltage  and  current,  the  power  and  energy  output  of  the  capillary  was  subsequently 
calculated.  These  data  are  shown  in  figure  3.  Typical  pressure-time  data  recorded  in  the 
simulator  are  shown  in  figure  4. 


3.  NGEN  Code  Description 


The  Army’s  NGEN3  code  is  a  multi-dimensional,  multi-phase  CFD  code  that  incorporates  three- 
dimensional  continuum  equations  along  with  auxiliary  relations  into  a  modular  code  structure 
(25-28).  Since  accurate  charge  modeling  involves  flow  field  components  of  both  a  continuous 
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Time  (ms) 


Figure  3.  Electrical  profiles  (current  in  amperes:  red;  voltage  in  volts:  blue; 

power  in  kilowatts:  green;  energy  in  joules:  brown)  used  in  the  ARL 
25 -mm  simulator  tests  (22). 


Figure  4.  Pressure-time  history  measured  in  the  ARL  2 5 -mm  simulator  for  (a)  0.4-mm  disk  gap  shown  on 
left  and  (b)  0.8-mm  disk  gap  shown  on  right  (21-23). 
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and  a  discrete  nature,  a  coupled  Eulerian-Lagrangian  approach  is  utilized.  On  a  sufficiently 
small  scale  of  resolution  in  both  space  and  time,  the  components  of  the  flow  are  represented  by 
the  balance  equations  for  a  multicomponent  reacting  mixture  describing  the  conservation  of 
mass,  momentum,  and  energy.  A  macroscopic  representation  of  the  flow  is  adopted  using  these 
equations  derived  by  a  formal  averaging  technique  applied  to  the  microscopic  flow.  These 
equations  require  a  number  of  constitutive  laws  for  closure  including  state  equations, 
intergranular  stresses,  and  interphase  transfer.  The  numerical  representation  of  these  equations 
as  well  as  the  numerical  solution  thereof  is  based  on  a  finite-volume  discretization  and  high- 
order  accurate,  conservative  numerical  solution  schemes.  The  spatial  values  of  the  dependent 
variables  at  each  time  step  are  determined  by  a  numerical  integration  method  denoted  the 
Continuum  Flow  Solver  (CFS),  which  treats  the  continuous  phase  and  certain  of  the  discrete 
phases  in  an  Eulerian  fashion.  The  Flux-Corrected  Transport  scheme  (29)  is  a  suitable  basis  for 
the  CFS  since  the  method  is  explicit  and  has  been  shown  to  adapt  easily  to  massively  parallel 
computer  systems.  The  discrete  phases  are  treated  by  a  Eagrangian  formulation,  denoted  the 
Large  Particle  Integrator  (LPI),  which  tracks  the  particles  explicitly  and  smoothes  discontinuities 
associated  with  boundaries  between  propellants  yielding  a  continuous  distribution  of  porosity 
over  the  entire  domain.  The  manner  of  coupling  between  the  CFS  and  the  LPI  is  through  the 
attribution  of  properties  (e.g.,  porosity,  and  mass  generation).  The  size  of  the  grid  as  well  as  the 
number  of  Lagrangian  particles  is  user  prescribed. 

For  the  simulations  of  novel  solid  propellant  configurations,  such  as  disks  stacked  axially  along 
the  chamber  centerline  and/or  thin  annular  concentric  layers  (wraps),  the  NGEN  code  takes  a 
macroscopic  approach.  These  solid  propellant  media  are  modeled  using  Lagrange  particles  that 
regress,  produce  combustion  product  gases,  and  respond  to  gas  dynamic  and  physical  forces. 
Individual  grains,  sticks,  slab,  and  wrap  layers  are  not  resolved;  rather,  each  medium  is 
distributed  within  a  specified  region  in  the  gun  chamber.  The  constitutive  laws  that  describe 
interphase  drag,  form- function,  etc.,  assigned  to  these  various  media,  determine  preferred  gas 
flow  paths  through  the  media  (e.g.,  radial  for  disks  and  axial  for  wraps)  and  responses  of  the 
media  to  forces.  Media  regions  can  be  encased  in  impermeable  boundaries  that  yield  to  gas 
dynamic  flow  after  a  prescribed  pressure  load  is  reached.  Details  of  the  NGEN  code  are  supplied 
elsewhere  (27,  28). 

For  modeling  ETC  igniters,  the  NGEN  code  is  linked  to  the  ARL  plasma  capillary  code 
developed  by  Powell  and  Zielinski  (30).  This  capillary  code  uses  the  length,  diameter,  and 
material  composition  of  the  ablation  capillary  as  well  as  the  current-time  profde  and  supplies  as 
output  the  time-dependent  temperature,  density,  molecular  weight,  and  flow  rate  of  the  plasma 
efflux.  This  information  is  mapped  to  the  computational  cells  that  lie  along  the  breech  face  of 
the  gun  chamber,  from  the  centerline  to  a  predetermined  radial  position.  Pressure  is  linked 
between  the  NGEN  code  and  the  capillary  code  so  that  the  changing  chamber  pressure,  as 
computed  by  the  NGEN  code,  is  input  the  capillary  code  and  modifies  the  plasma  efflux.  The 
NGEN  code  represents  the  plasma  as  a  single  component  gas  with  properties  supplied  by  the 
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capillary  code.  Of  course,  the  plasma  is  actually  a  chemically  diverse  and  reactive  media  which 
has  been  treated  as  such  in  separate  modeling  studies  (17-19).  The  macroscopic  approach  taken 
in  the  NGEN  code,  with  respect  to  plasma  representation,  has  been  used  to  effectively  model 
several  solid  propellant  charge  configurations  (24,  31-32). 


4.  Modeling  Results  for  the  ARL  25-mm  Simulator 


4.1  Simulations  for  Inert  Disks 

In  a  previous  paper  (24),  the  NGEN  code  was  used  to  simulate  the  case  of  inert  (Masonite)  disks 
placed  in  the  25-mm  simulator  in  which  high  temperature  plasma  was  injected.  This  setup  is 
described  in  the  experimental  description  section  of  this  report  and  by  Nusca  and  Howard  (24). 
Comparison  between  measured  and  computed  pressure  data  for  that  work  demonstrated  the 
ability  of  the  NGEN  code  to  capture  the  plasma  injection  and  dispersal  processes  within  a 
chamber  tightly  loaded  with  propellant  disks  but  did  not  test  the  ability  of  the  code  to  simulate 
propellant  ignition  and  flame  spreading.  Figure  5  shows  the  computed  results  (24)  for  the  case 
of  0.4-  and  0.8-mm  gaps  between  the  inert  disks.  These  results  can  be  compared  to  experimental 
data  in  figure  4  (note  that  a  common  color  code  has  been  used  so  that  data  from  the  pressure  tap 
located  near  the  breech  end  on  the  chamber  are  colored  black  and  data  from  the  pressure  tap 
located  in  the  chambrage  region,  loaded  with  the  inert  grains,  are  colored  red).  The  NGEN 
simulations  feature  a  pressure  data  station  midway  between  the  rear  and  the  forward  stations, 
colored  green,  that  was  not  used  in  the  simulator.  Both  the  measured  and  computed  pressure  for 
the  rear  location  show  an  initial  rise  starting  before  0.1  ms  and  peaking  at  about  0.35  ms  with  a 
pressure  of  about  1 . 1  MPa.  The  pressure  traces  for  the  forward  location  show  a  similar  behavior 
with  a  peak  of  about  0.9  MPa.  As  is  expected,  each  pressure  tap  included  in  the  simulation  is 
affected  in  turn,  indicative  of  a  pressure  wave  that  travels  through  the  chamber  from  breech  end 


Figure  5.  Computed  pressure-time  history  for  (a)  0.4-mm  disk  gap  shown  on  the  left  and  (b)  0.8-mm  disk 
gap  shown  on  the  right  (24). 
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to  forward  chambrage;  recall  from  figure  3  that  the  peak  current  occurs  at  about  0.2  ms  and  that 
the  capillary  is  no  longer  charged  after  about  0.38  ms.  While  the  measured  data  traces  show  a 
drop  in  pressure  to  a  level  of  about  0.6  MPa  after  about  0.35  ms,  the  computations  show  that  the 
pressure  in  the  chamber,  for  all  tap  locations,  reaches  a  steady  pressure  of  about  1  MPa.  It  is 
thought  that  some  loss  mechanism  particular  to  this  experiment  is  missing  from  the  NGEN 
model  or  that  the  simulator  was  prone  to  developing  leaks.  To  date  these  issues  have  not  been 
specifically  resolved,  but  it  was  thought  that  the  comparisons  for  live  propellant  disks  (JA2) 
should  proceed  in  the  current  report. 

4.2  Simulations  for  JA2  Disks 

The  experimental  setup  for  the  case  of  JA2  propellant  disks  loaded  into  the  same  simulator  used 
for  inert  disk  experiments  is  described  in  the  experimental  description  section  of  this  report.  The 
NGEN  code  was  setup  as  closely  as  possible  to  resemble  the  experiment.  Figure  6  shows  a 
series  of  four  pulses  delivered  to  the  polyethylene  capillary  in  order  to  produce  four  distinct 
plasma  jets  that  issue  into  the  propellant  filled  chamber.  For  a  “low  energy”  of  about  500  J 
(actually  260-680  J)  both  a  short  and  a  long  pulse  of  about  0.3  and  1  ms,  respectively,  were  used. 
For  a  “high  energy”  of  about  1000  J  (actually  780-1470  J)  both  a  short  and  a  long  pulse  of  about 
0.3  and  1  ms,  respectively,  were  also  used.  These  data  are  used  in  the  NGEN  code  to  form  the 
plasma  efflux  into  the  chamber,  as  previously  discussed. 


Figure  6.  Experimental  current  profiles  derived  from  the  ARL  25-mm 
simulator  tests  and  used  in  the  NGEN  code  simulations. 


7 


Figure  7  shows  the  computed  pressure  results  for  the  rear  position  (i.e.,  location  of  pressure  tap 
placed  along  the  chamber’s  radial  wall  and  near  the  rear  wall)  and  the  forward  position  (i.e., 
location  of  pressure  tap  placed  along  the  chamber’s  radial  wall  and  into  the  forward  chambrage) 
using  the  short  duration  plasma  pulse  with  both  low  and  high  energies  (values  as  previously 
described).  Immediately  evident  is  that  the  pressure  never  rises  above  about  3.5  MPa  for  the  low 
energy  pulse  (figure  7a)  but  exceeds  the  yield  pressure  of  the  simulator  (about  13  MPa)  for  the 
high  energy  pulse  (figure  7b).  Thus,  it  can  be  said  that  both  the  experiments  and  the  simulations 
indicate  that  the  low  energy  pulse  is  not  sufficient  to  ignite  the  JA2  propellant  disks,  whereas  the 
high  energy  pulse  is  more  than  sufficient  to  induce  propellant  ignition.  For  the  simulations,  the 
chamber  pressure  continues  to  rise  slowly  for  the  case  of  no  propellant  ignition  because  heat  loss 
of  the  chamber  wall  is  not  included  in  the  simulation  and  rises  well  above  20  MPa  for  the  case  of 
propellant  ignition  since  chamber  yield  is  also  not  included  in  the  simulation.  Both  the  computed 
and  measured  pressures  for  the  rear  position  rise  rapidly,  followed  by  pressurization  of  the 
forward  end  of  the  chamber  after  some  time  delay.  Both  the  computed  and  measured  pressure 
traces  show  significant  pressure  wave  action  and  oscillations  at  early  times.  The  simulations 
seem  to  overprotect  both  the  early-time  peak  pressure  at  the  rear  position  and  the  rate  at  which 
the  pressure  wave  reaches  the  forward  end  of  the  chamber.  In  addition,  the  simulations  predict 
an  earlier  than  measured  propellant  ignition  (i.e.,  pressure  runaway)  for  the  high  energy  pulse 
(figure  7b). 


Short  Pulse  with  High  Energy 

COMPUTED  -  BEAR  POSITION 


Figure  7.  Pressure  tap  results  for  JA2  propellant  disks  in  the  25-mm  simulator  using  a  short  time 

duration  current  pulse  (a)  for  the  low  energy  pulse  as  shown  on  the  left  and  (b)  for  the  high 
energy  pulse  shown  on  the  right. 

Figure  8  shows  the  computed  pressure  results  for  the  rear  position  and  the  forward  position  using 
the  long  duration  plasma  pulse  with  both  low  and  high  energies  (values  as  previously  described). 
Immediately  evident  is  that  the  pressure  never  rises  above  about  4.5  MPa  for  the  low  energy 
pulse  (figure  8a)  but  exceeds  the  yield  pressure  of  the  simulator  (about  13  MPa)  for  the  high 
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Figure  8.  Pressure  tap  results  for  JA2  propellant  disks  in  the  25 -mm  simulator  using  a  long  time  duration 

current  pulse  (a)  for  the  low  energy  pulse  as  shown  on  the  left  and  (b)  for  the  high  energy  pulse  shown 
on  the  right. 

energy  pulse  (figure  8b).  Thus,  it  can  be  said  that  both  the  experiments  and  the  simulations 
indicate  that  the  low  energy  pulse  is  not  sufficient  to  ignite  the  JA2  propellant  disks,  whereas  the 
high  energy  pulse  is  more  than  sufficient  to  induce  propellant  ignition.  Both  the  computed  and 
measured  pressures  for  the  rear  position  rise  rapidly,  followed  by  pressurization  of  the  forward 
end  of  the  chamber  after  some  time  delay.  Both  the  computed  and  measured  pressure  traces 
show  significant  pressure  wave  action  and  oscillations  at  early  times.  The  simulations  seem  to 
overprotect  both  the  early-time  peak  pressure  at  the  rear  position  and  the  rate  at  which  the 
pressure  wave  reaches  the  forward  end  of  the  chamber.  The  simulation  predicts  a  gradual 
transition  to  propellant  ignition  for  the  high  energy  pulse  (figure  8b)  in  close  agreement  with  the 
measurements. 

The  close  agreement  between  measured  and  computed  pressures  for  the  longer  and  more 
energetic  plasma  events  would  seem  to  indicate  that  physical  models  of  the  plasma  properties 
(e.g.,  thermodynamic  equations  and  the  equation  of  state  [EOS])  resident  in  the  NGEN  code 
require  improvement.  When  the  plasma  is  more  slowly  injected  into  the  chamber  (i.e.,  long 
duration  pulse)  and  propellant  ignition  occurs  the  gases  in  the  chamber  are  more  closely 
simulated  using  a  mixture  of  high-temperature  air  and  reacted  JA2  with  relatively  little  plasma. 
Since  the  EOS  for  the  plasma  is  probably  the  least  accurate,  the  result  is  a  more  accurate  overall 
simulation.  Another  explanation  of  the  difference  in  behavior  of  the  model  with  the  results  in 
figures  7b  and  8b  is  the  amount  of  energy  deposited  in  the  JA2-filled  chamber  as  a  function  of 
pulse  length.  Figure  9  shows  the  maximum  pressure  for  each  experimental  test  plotted  as  a 
function  of  plasma  energy  for  that  particular  test.  For  the  experimental  tests  where  insignificant 
JA2  combustion  occurred,  the  pressures  should  be  a  function  of  the  energy  utilized  by  the  plasma 
capillary.  When  these  pressures  are  plotted  as  in  figure  9  this  function  appears  to  be  linear. 

When  the  pressures  in  figures  7b  and  8b  that  occur  just  prior  to  the  onset  of  ignition  are  included 
in  figure  9,  the  pressure  generated  by  the  short  pulse  appears  to  be  below  that  predicted  by  the 
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♦  Short  Pulse 
■  Long  Pulse 


Energy  (J) 


Figure  9.  Plot  of  maximum  pre-ignition  pressure  for  25-mm  simulator  filled 
with  JA2  disks  as  a  function  of  plasma  energy  in  joules. 

linear  function  and  probably  indicates  that  the  ignition  in  figure  7b  is  weaker  than  expected.  If 
the  same  pressure  is  plotted  for  the  long  pulse,  it  is  higher  than  that  predicted  by  the  linear 
assumption.  Therefore,  it  appears  that  more  than  sufficient  energy  is  present  for  ignition  in  this 
case.  The  model  and  experiment  should  be  somewhat  in  agreement  as  shown  in  figure  8b.  If,  in 
the  case  of  the  short  pulse,  the  trend  to  ignition  is  weaker  than  the  model  is  expecting,  the  run  up 
to  ignition  in  the  experiment  would  be  delayed  relative  to  that  of  the  model. 

A  more  detailed  examination  of  propellant  pressurization  and  propellant  ignition  follows;  the 
reader  is  directed  to  Nusca  and  Howard  (24)  for  a  discussion  of  NGEN  code  simulations  for  the 
case  of  no  propellant  ignition,  as  is  the  case  for  the  low  energy  pulse  in  the  current  work. 

Figure  10  shows  a  sequence  of  eight  color  flowfield  pressure  maps  computed  as  part  of  the 
NGEN  simulation  from  0.22  to  2.5  ms  for  the  case  of  a  long  duration  current  pulse  of  high 
energy,  corresponding  to  figure  8b.  A  constant  pressure  range  (i.e.,  blue  is  0.1  MPa  and  below 
while  red  is  4  MPa  and  above)  was  used  for  the  plots  from  0.22  to  0.90  ms  and  then  a  larger 
pressure  range  (i.e.,  blue  is  0.1  MPa  and  below  while  red  is  12  MPa  and  above)  was  used  for  the 
plots  from  0.98  to  2.5  ms.  This  was  done  because  after  about  1  ms  the  chamber  is  essentially  of 
constant  pressure  that  is  steadily  increasing  (recall  figure  8b).  Superimposed  on  the  pressure 
fields  are  velocity  vectors  shown  in  black  that  indicate  the  magnitude  and  direction  of  the  gas 
flow  (note  that  not  all  vectors  are  plotted  for  clarity).  The  velocity  vectors  illustrate  that  the 
plasma  efflux  is  emanating  from  a  small  region  along  the  breech  wall  and  near  the  chamber 
centerline.  It  can  be  seen  that  both  the  empty  center  core  of  the  chamber  and  the  narrow  space 
along  the  breech  wall,  that  are  not  occupied  by  propellant,  pressurize  rapidly,  followed  by  a 
pressure  wave  that  moves  from  rear  to  forward  through  the  disks  (note  the  velocity  vectors). 


10 


Figure  10.  Color  pressure  contours  and  superimposed  velocity  vectors  as  computed  using 
the  NGEN  code  for  the  case  of  JA2  disks  in  a  2 5 -mm  simulator,  long  duration 
pulse  with  high  energy,  and  for  times  from  0.22  to  2.5  ms. 
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When  this  wave  reaches  the  base  of  the  projectile  afterbody  (-0.68  ms),  the  chambrage  region  is 
pressurized,  the  wave  reflects  and  then  travels  rearward  toward  the  breech.  By  about  1  ms,  the 
plasma  efflux  has  terminated,  most  pressure  wave  action  has  ceased,  and  the  chamber  is 
pressurized  to  near  uniformity.  Subsequently,  the  chamber  pressure  continues  to  climb  because 
the  propellant  disks  have  ignited  and  are  generating  combustion  gases. 

Figure  1 1  shows  a  sequence  of  eight  color  propellant  temperature  maps  computed  as  part  of  the 
NGEN  simulation  from  0.22  to  2.5  ms  for  the  case  of  a  long  duration  current  pulse  of  high 
energy,  corresponding  to  figures  8b  and  10.  In  this  case,  the  color  blue  represents  a  region  that  is 
either  not  occupied  with  propellant  or  a  region  that  is  occupied  with  propellant  at  294  K.  The 
color  red  indicates  a  region  of  fully  ignited  propellant  while  colors  between  blue  and  red  indicate 
propellant  that  is  being  heated  but  has  yet  to  ignite.  Superimposed  on  these  fields  are  velocity 
vectors  shown  in  black  that  indicate  the  magnitude  and  direction  of  the  gas  flow  (note  that  not  all 
vectors  are  plotted  for  clarity).  The  velocity  vectors  illustrate  that  the  plasma  efflux  is  emanating 
from  a  small  region  along  the  breech  wall  and  near  the  chamber  centerline  that  rapidly  heats  and 
ignites  the  edge  of  the  propellant  disks  immediately  adjacent  to  this  efflux.  As  the  pressure  wave 
moves  forward  through  the  empty  center  core  of  the  chamber  and  plasma  is  forced  between  the 
disks  (note  vectors  for  0.84  ms),  propellant  ignition  is  started.  Note  that  the  inert  propellant 
grains  in  the  forward  chambrage  region  are  heated  to  ignition  temperature  (0.84  ms)  but  they  are 
not  allowed  to  either  ignite  or  generate  combustion  gases  as  part  of  the  simulation.  As  the 
pressure  wave  moves  rearward  (after  0.68  ms)  the  propellant  disks  are  slowly  ignited  first  along 
the  edges  adjacent  to  the  chamber  center  core  and  then  through  the  propellant  as  convective  heat 
transfer  is  accomplished  by  upward  moving  plasma.  Most  wave  action  in  the  chamber  ceases 
after  0.98  ms  and  the  entire  bed  of  propellant  disks  are  fully  ignited  by  2.5  ms.  These 
observations  correspond  to  the  pressure  traces  of  figure  8b  and  explain  the  smooth  and  gradual 
rise  in  pressure  after  a  steady-state  is  reached  (i.e.,  pressure  waves  have  ceased)  and  the 
propellant  is  ignited  and  generating  propellant  gas. 

Figures  12  and  13  show  a  sequence  of  eight  color  flowfield  pressure  maps  and  eight  color 
propellant  temperature  maps,  respectively,  computed  as  part  of  the  NGEN  simulation  from 
0.22  to  2.5  ms  for  the  case  of  a  long  duration  current  pulse  of  low  energy  (see  figure  8a).  A 
constant  pressure  range  (i.e.,  blue  is  0.1  MPa  and  below  while  red  is  4  MPa  and  above)  was  used 
in  figure  12  for  plots  from  0.22  to  0.90  ms  and  then  a  larger  pressure  range  (i.e.,  blue  is  0.1  MPa 
and  below  while  red  is  12  MPa  and  above)  was  used  for  the  plots  from  0.98  to  2.5  ms.  For 
figure  13,  the  color  blue  represents  a  region  that  is  either  not  occupied  with  propellant  or  a  region 
that  is  occupied  with  propellant  at  294  K,  while  the  color  red  indicates  a  region  of  fully  ignited 
propellant;  colors  between  blue  and  red  indicate  propellant  that  is  being  heated  but  has  yet  to 
ignite.  Superimposed  on  the  pressure  and  temperature  fields  are  velocity  vectors  shown  in  black 
that  indicate  the  magnitude  and  direction  of  the  gas  flow  (note  that  not  all  vectors  are  plotted  for 
clarity). 
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PROPELLANT  TEMP.  CONTOURS  (294  <  T  <  444  K)  &  SELECTED  VEL.  VECTORS 
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Figure  11.  Color  propellant  temperature  contours  and  superimposed  velocity  vectors  as 

computed  using  the  NGEN  code  for  case  of  JA2  disks  in  simulator,  long  duration 
pulse  with  high  energy,  and  for  times  from  0.22  to  2.5  ms. 
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Figure  12.  Color  pressure  contours  and  superimposed  velocity  vectors  as  computed  using 
the  NGEN  code  for  the  case  of  JA2  disks  in  a  2 5 -mm  simulator,  long  duration 
pulse  with  low  energy,  and  for  times  from  0.22  to  2.5  ms. 
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Figure  13.  Color  propellant  temperature  contours  and  superimposed  velocity  vectors  as 
computed  using  the  NGEN  code  for  case  of  JA2  disks  in  simulator,  long 
duration  pulse  with  low  energy,  and  for  times  from  0.22  to  2.5  ms. 
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Comparing  the  results  of  the  high  energy  pulse  (figures  10  and  1 1)  to  those  for  the  low  energy 
pulse  (figures  12  and  13),  it  can  be  observed  that  after  the  plasma  efflux  ceases  (i.e.,  after  about 
0.98  ms)  a  high  energy  pulse  causes  propellant  ignition  while  the  low  energy  pulse  does  not.  For 
the  low  energy  pulse,  pressure  levels  between  the  disk  propellant,  especially  in  the  forward 
section  of  the  chamber  are  too  low  to  create  enough  convective  heat  transfer  into  the  propellant 
for  ignition,  suggesting  that  a  link  exists  between  the  energy  content  of  the  plasma  and  propellant 
ignition  which  is  similar  to  the  case  of  conventional  hot  gas  igniters.  It  is  clear  that  the  presence 
of  plasma,  per  se,  does  not  guarantee  propellant  ignition,  but  that  the  plasma  efflux  must  be 
sufficiently  energetic  regardless  of  the  duration. 

These  points  are  reinforced  by  the  results  displayed  in  figures  14  and  15,  showing  a  sequence  of 
eight  color  propellant  temperature  maps  computed  as  part  of  the  NGEN  simulation  from 
0.22  to  2.5  ms  for  the  case  of  a  short  duration  current  pulse  of  either  high  or  low  energy, 
respectively  (recall  figure  7).  The  high  energy  pulse  of  short  duration  rapidly  ignites  the 
propellant  charge  (even  more  rapidly  than  the  same  pulse  but  of  long  duration  [figure  11])  while 
the  low  energy  pulse  or  short  duration  merely  heats  the  edge  of  the  propellant  disks  and  fails  to 
effectively  penetrate  the  charge. 


5.  Summary  and  Conclusions 


The  NGEN  multidimensional,  multiphase  CFD  code  has  been  successfully  used  to  model  the 
interior  ballistics  of  solid  propellant  charges  ignited  using  hot  chemical  gases  like  those 
generated  by  burning  black  powder  or  Benite  (27,  28).  In  order  to  model  charges  ignited  using 
injected  plasma,  which  is  generated  by  an  ablation  capillary  of  polyethylene,  the  NGEN  code 
was  upgraded  with  a  single-component,  non-reacting  gas  phase  treatment  of  the  plasma  and 
linked  to  a  well-known  and  independently  validated  capillary  model. 

In  previous  works,  NGEN  simulations  for  plasma  injection  into  the  gun  chamber  showed  that  the 
ignition  delay  is  lengthened,  even  to  the  point  of  a  non- ignition  event,  when  the  current  that  is 
input  to  the  ablation  capillary  is  reduced  (24,  31,  32).  This  is  an  important  result  as  it 
necessitates  a  closer  linkage  between  the  charge  designer  and  the  igniter  designer  when  the  ETC 
option  is  utilized.  Given  independent  plasma-propellant  modeling  results  that  point  to  the  need 
for  a  multicomponent,  reacting  flow  representation  of  the  plasma  (17-19)  and  the  incompatibility 
of  this  level  of  plasma  modeling  with  the  NGEN  code,  within  the  constraints  of  currently 
available  computers,  there  is  an  immediate  need  to  validate  the  NGEN  code  for  the  design  of 
ETC  guns. 
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PROPELLANT  TEMP.  CONTOURS  (294  <  T  <  444  K)  &  SELECTED  VEL.  VECTORS 
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Figure  14.  Color  propellant  temperature  contours  and  superimposed  velocity  vectors  as 
computed  using  the  NGEN  code  for  case  of  JA2  disks  in  simulator,  short 
duration  pulse  with  high  energy,  and  for  times  from  0.22  to  2.5  ms. 
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PROPELLANT  TEMP.  CONTOURS  (294  <  T  <  444  K)  &  SELECTED  VEL.  VECTORS 


Figure  15.  Color  propellant  temperature  contours  and  superimposed  velocity  vectors  as 
computed  using  the  NGEN  code  for  case  of  JA2  disks  in  simulator,  short 
duration  pulse  with  low  energy,  and  for  times  from  0.22  to  2.5  ms. 
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A  significant  step  in  this  effort  has  been  to  compare  the  chamber  pressure  measured  in  a  sub¬ 
scale  ballistic  simulator  with  that  computed  using  the  code.  The  configuration  chosen  was  a 
25 -mm  chamber  filled  with  either  inert  disks  or  disks  made  from  JA2  propellant  and  linked  to  an 
ablative  capillary  of  the  size  and  power  of  a  typical  ETC  gun  igniter.  Previous  works  showed 
that  the  NGEN  code  was  able  to  capture  the  essential  features  of  the  plasma  efflux  as  it  was 
injected  and  then  spread  around  the  inert  disks  via  diffusion/convection  (24).  In  the  present 
work,  using  a  chamber  filled  with  JA2  disks,  the  NGEN  code  results  showed  that  a  low  energy 
pulse  (500  J),  whether  short  (~0.3  ms)  or  long  (~1  ms)  in  duration,  was  not  sufficient  to  ignite  the 
disks,  while  a  high  energy  pulse  (1000  J)  of  either  duration  was  able  to  accomplish  JA2  ignition. 
These  results  are  in  full  agreement  with  the  experiment.  Detailed  comparison  between  the 
measured  and  computed  pressure  histories  showed  some  discrepancies  that  have  been  attributed 
to  the  state  equations  employed  by  the  model  and/or  an  apparent  ignition  “threshold”  in  terms  of 
plasma  energy  that  may  have  biased  the  comparison.  A  remarkable  agreement  between 
measured  and  computed  pressure  history  was  achieved  for  the  case  of  a  long  duration  pulse  and 
high  energy. 
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1  CDR  USAATC 
CSTE  DTC  AT  SL 
APG  MD  21005 

66  DIR  USARL 

AMSRD  ARL  WM 
T  ROSENBERGER 
AMSRD  ARL  WM  M 
S  MCKNIGHT 
AMSRD  ARL  WM  T 
B  BURNS 

AMSRD  ARL  WM  TB 
P  BAKER 

AMSRD  ARL  WM  B 
C  CANDLAND 
J  MORRIS 
J  NEWILL 
M  ZOLTOSKI 
AMSRD  ARL  WM  BA 
B  DAVIS 
G  BROWN 
D  HEPNER 
G  KATULKA 
T  KOGLER 
D  LYON 

AMSRD  ARL  WM  BC 
M BUNDY 
G  COOPER 
J  DESPIRITO 
J  GARNER 
P  PLOSTINS 
JSAHU 
S  SILTON 
P  WEINACHT 
AMSRD  ARL  WM  BD 
W  ANDERSON 
R  BEYER 
A  BRANT 
SBUNTE 
L  CHANG 
T  COFFEE 
J  COLBURN 
P  CONROY 
N  ELDREDGE 
B  FORCH 
B  HOMAN 
A  HORST 

S  HOWARD  (6  CPS) 
PKASTE 
A  KOTLAR 
C  LEVERITT 
RLIEB 


K  MCNESBY 
M  MCQUAID 
A  MIZIOLEK 
J  NEWBERRY 
M  NUSCA  (6  CPS) 

R  PESCE-RODRIGUEZ 
S  PIRIANO 
G  REEVES 
BRICE 
R  SAUSA 
J  SCHMIDT 
A  WILLIAMS 
AMSRD  ARL  WM  BF 
R  ANDERSON 
W  OBERLE 
D  WILKERSON 
AMSRD  ARL  WM  EG 
E  SCHMIDT 
AMSRD  ARL  WM  SG 
W  CIEPIELA 
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Intentionally  left  blank. 
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